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We consider a SU (N) x SU{M) generalization of the multichannel single-impurity 
Kondo model which we solve analytically in the limit N oo, M — > oo, with 
7 = M/N fixed. Non-Fermi liquid behavior of the single electron Green function 
and of the local spin and flavor susceptibilities occurs in both regimes, N < M and 
N > M, with leading critical exponents identical to those found in the conformal 
field theory solution for all N and M (with M > 2). We explain this remarkable 
agreement and connect it to "spin-flavor separation", the essential feature of the 
non-Fermi-liquid fixed point of the multichannel Kondo problem. 
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In recent years, the multichannel Kondo model, first introduced by Nozieres and 
Blandin has been the focus of intense activity. This model involves a single local mo- 
ment of spin Si antiferromagnetically coupled to M identical conduction bands in an A^<j-site 
system, as described by the Hamiltonian, 

H = 5Z (^kcl^aCkaa + {J/Ns)Sl ■ [ XI 4 

kaa ka,k'cr' Q? 

where c\^^ creates a conduction electron of (radial) wave vector k, spin cr =|, | and channel 
index a running from 1,..,M; the coupling > is antiferromagnetic. It is now estab- 
lished that in the "overcompensated" regime, 25*/ < M, the screening of the local 
moment by the conduction electrons drives the metal into a non-Fermi liquid critical state 
at T = 0, fully validating Nozieres and Blandin's original arguments. 

The renewed interest in this model is, in part, due to a number of suggestions for real- 
izations of two-channel overcompensated behavior in explicit experimental contexts, most 
notably: (i) the quadrupolar Kondo effect ||^ in heavy fermion alloys (ii) non- Fermi 
liquid scattering rates in narrow copper point-contacts (iii) "marginal" Fermi liquid nor- 
mal state properties ^ of the high-Tc materials Moreover, it has been argued that 
the two-channel Kondo model provides a link to exotic superconductivity Related 
non-Fermi liquid behavior and enhanced pairing correlations have been recently found in 
the mixed-valence regime at low temperatures in certain extended Anderson single-impurity 
models motivated by the electronic structure of the high-Tc materials . 

From a theoretical point of view, remarkable progress has been made in understanding 
the universal properties of the multi-channel Kondo and other impurity models through 
the use of conformal field theory techniques |^. This approach provides a clear picture 
of the separation of spin, channel, and charge excitations at non-Fermi liquid fixed points 
and describes the subtle recombination of these degrees of freedom required in the Fermi 
liquid case. In spite of its elegance and power this method only classifies the possible critical 
behaviors without providing a constructive route to the solution of a particular model. 

In this letter we formulate a controlled calculation method which can, in principle, in- 
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corporate both the comphcations of real materials and the conceptual insights gained from 
conformal field theory. Apart from the obvious practical application to dilute impurity 
systems, this question is important in studying the possibility for non-Fermi liquid behav- 
ior in lattice systems, for which the conformal field theory techniques are not immediately 
applicable. An explicit route for addressing the latter problem, which appears particularly 
promising, involves constructing mean-field theories for the lattice by solving single-impurity 
models embedded in a self-consistent medium |llO[| . 

Below we concentrate on the SU{N) x SU (M) generalization of the multichannel Kondo 
Hamiltonian of Eq. (1), where and M are the degeneracies in the spin and flavor quantum 
numbers, respectively. We use a functional integral approach based on the "slave-Boson" 
representation [0] which explicitly separates the (local) spin and flavor excitations. The 
limit A^, M — > oo with 7 = M/N fixed leads to a closed set of coupled self-consistent integral 
equations which can be solved analytically in the asymptotic low-frequency, zero temper- 
ature limit. These are identical to the "Non-Crossing Approximation" (NCA) equations of 
perturbation theory |TT|. Although these equations break down for sufficiently low energies 



in the Fermi liquid (single channel) case []T5|, it is the unique feature of our large N,M 



treatment of the multi-channel Kondo problem that the NCA becomes exact. 

In fact, we show that the single electron Green function, and the local spin and flavor 
susceptibilities display non-Fermi liquid behavior with leading critical exponents identical 
to those found in the conformal field theory solution for all N,M > 2. We explain this 
apparently surprising result by demonstrating that the fiuctuations provide no corrections 
to the leading exponents obtained from NCA to all orders inl/N,l/M |T^. Our calculation 
is an explicit realization of the "spin-charge separation" idea emphasized by Anderson in his 
theory of the high-Tc materials based on the single-band Hubbard model [0. Ultimately, 
we hope that the considerations presented here will be useful in treating lattice Fermion 
systems with non-Fermi liquid ground states. 

Our starting point is a path integral treatment of the generalized infinite U Anderson 
model Hamiltonian, 



H=Y: e,4,,,„c,,.,„ + efY: flu + {V/^N,) ^ [f%c,^^^^ + h.c] (2) 

where the Fermion, /J, creates a local spin excitation and the Boson, 6^, transforms according 
to the conjugate representation of SU{M), and annihilates the flavor quantum number of 
the "vacuum" state produced by destroying a conduction electron. The "completeness" of 
the local states at the impurity site, represented by the constraint, J2a flfa + J2a baba = 1, 
is implemented in the usual way [jll|, by introducing a fictitious field A coupling to the 
constrained charges, which is taken to oo at the end of the calculation. In the limit e/ < 
0,V/\ef \ < 1 (D leads to the SU{N) x SU{M) Coqblin-Schrieffer model, with exchange 
coupling J = V'^/\ef\. 

The most compact way of presenting our arguments is in terms of the impurity contri- 
bution to the partition function, Zimp, obtained after performing a Gaussian integral over 
the conduction electron fields. (The latter contribute an overall multiplicative factor, Zc, 
of the free-electron gas partition function to the full partition function, Z = ZcZi„ip.) By 
completing squares, the hybridization term is eliminated in favor of a spin-flavor interac- 
tion contribution, S^^^t = -iVyN)^^,^ J dr J dr' j2{T)f,{r')G'{r - T')bUr')b^ir), where 
G^{t — t') = —J2k{d/dT + ek)^^/Ns is the non-interacting conduction electron Green 
function at the impurity site; and, as usual, in order to obtain a nontrivial large N, M 
limit, we have defined a rescaled hybridization matrix element, V = \/NV, which should 
be considered of order unity at the end of the calculation. The next step is to intro- 
duce two composite fields, non-local in imaginary time, ^f^{r,T') and ^^^(t, t') (with 
$j;,(r, r') = $/^5(r',r)), which decouple the interaction term, Smt, enabling us to write, 
Zimp = I[^f][^b][^^f][1^^b]GW{~S)- The effective action, 5*, is then given by 

^ = E dr'flir) [5{t - r')(|: + + A) + ^ E $,,(r, r')G°(r - r')]/.(r') 

-^Hfdr f dr'^j^{T',T)G\r-r')^Ur,r'). 
N ^Jo Jo 

Finally, the large N, M calculation proceeds as usual by carrying out the Gaussian integral 



over the local Fermions and Bosons, and ba, to produce an effective action for the com- 
posite fields, ^fa,ba', this is then evaluated by a saddle point integration (Hereafter we 
drop the bar on the channel index, a.) 

In the N, M = oo limit with 7 = M/N fixed the saddle point approximation be- 
comes exact and leads to time translationally invariant solutions, $/,b(T — t')- After 
Fourier transforming in terms of Fermionic and Bosonic imaginary frequencies, Un, i^n, 
the saddle point equations can be written as $/(ito'„) = [iun — ~ ^/(^^n)]"^ and 
^biit^n) = [i^n — nf,(ii/„)]^^ This defines the self energies, S/ and Ilf,, which, on the real 
frequency axis {iun = to + iO^, ivn = oj + ^0^), satisfy the self-consistent equations, 

T.f{uo) = ^— I dtf{t)H{t + uj) (4) 
vr J 

n,M = - / def{e)^f{e + u). (5) 

Here, F = TrpV"^ is the bare hybridization width (which includes a factor of introduced by 
the rescaling), and p is the conduction electron density of states, assumed constant in the 
energy range of interest. In principle, the functions Sj and H;, elsewhere in the complex plane 
can be obtained by the appropriate analytic continuation. Also note that in Eq. (^ the 
A 00 limit has already been taken, and that, due to the invariance under SU [N) x SU (M) 
rotations the spin and flavor indices have dropped out. 

In the low-frequency, low-temperature limit. Equations (|^) can be easily solved either 
by direct substitution or by reducing them to differential equations as done for the single 
channel case by Miiller-Hartmann [O. The relevant analysis can be found in a number of 



places in the literature ||Tl],|12[, and will not be repeated here. 

The saddle-point solutions for $ f and can be written in terms of the reduced frequency 
variable, = [((1 + 'j)/'j){Eq — cu) /Tq]^/^^^'''\ where Eq is the ground state energy and 
To = Z}(7r/7rD)'''exp(7rej/F) is the Kondo scale. The calculation is most easily done for 
negative frequencies, u < Eq, and the appropriate analytic continuation extends the results 
to > 0. At T — s> the resulting spin Fermion and flavor Boson spectral functions, 
Af^h{uj) = lm^f h{uj — zO^)/7r = A^j^fj {uj)6{uj — Eq) vanish for uj < Eq. In the process of 
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calculating physical quantities we will also need the spectral function for occupied states 
{uj < Eq), defined as A^j^j^{uj) = \imT^o[lni^ f^b{^ — iO^) / tt] exp{f3[Eo — uj]). Close to the 
threshold at Eq, these quantities take the form 



Afiuj) = sin(-^)|0|-^[l + 4-^ cos(-^)|e| + ...] (6) 
nTo 1 + 7 2 + 7 I+7 

1 • / ^7 M^i-ln . ^ch , T^l 



- ^ -(T^)|e|-[i - 4^ cos(^),er . ...] (7) 



= -Tfrioni - ^TT^m + •••] (8) 

■' iriQ 2 + 7 



4"^M = -|0r^[l -4^^161- + ...]. (9) 

Here, Wch = t^Tq/T is the weight of channel fluctuations in the ground state, and the 
constant, Z is related to the degeneracy of the impurity ground state per local degree of 
freedom (i.e., divided by A^(l + 7)) |T2|. Also, note the explicit breaking of particle-hole 
symmetry displayed by the positive (^/^'') and negative {A^fJ) contributions, consistent with 
the non-symmetric form of @). 

The scaling dimensions of the spin and flavor fields. A/ = [2(1 + 7)]^^ and A;, = 7A/ 
can be read off from the frequency dependence, \Eq — cjp^-f *"^, in Eq. (0). It then follows 
that the spin, channel and physical Fermion fields, all of which are bilinears of / and 6, have 
scaling dimensions, A^p = 2A/ = 1/(1+7), ^ch = 2Afc = 7/(1+7) and Ap = Af + Ab = 1/2, 
respectively. The resulting leading frequency dependence of the corresponding correlation 
function, \Eq — cup^sp.ch.F-i jg indeed what we observe. Below we summarize some of our 
explicit results. 

Single Electron Green Function: The local electron Green function, Qa,a, can be calcu- 
lated as a convolution of the spin and flavor propagators, leading to a local spectral function, 
Pa,ij.{uj,0) = Im ^o-,a (^^ — "^0^) /tt, of the form, 

PaA^,0)^ 7r/[(l + 7)' iV f] [l + 9iu)UiCu)+9i-co)f4Cu)], (10) 

with f±{u) = (a±p|^=f + 6±|cD|^-''), a_ = - [47/ (2 + 7) vr] sin(7rAeh)5(2A,p, A^,,), a+ = 
-cos(7rAc/,)a_, 6+ = -[4PVc/i/(l + 27)7?] sin(7rAc/i)5(2Ac/i, A^p) and 6_ = cos(7rAc/i)6+. 



Here a) = [(1 + 7)/7](ct;/To) and B{x, y) is the Beta function. Note that, in the overscreened 
case, M > N, the leading frequency dependence is the same as that obtained by Affleck and 



Resistivity: From the single-particle Green function we can obtain the resistivity from 
the standard transport theory formula, p{T) ~ [/ (ie(^^)r(e, T)]"-*^, in terms of the 
scattering rate, Ta^(uj,T)^^ = —2lmt'^^^{uj + iO'^,T) = [2Tpa^{uj,T)]/{pN), where we 
have used the relation between the conduction electron t-matrix and the single-electron 
Green function, t^^{LU,T) = V'^Qa^ii'^iT). The resulting leading behavior, p{T)/p{0) ~ 
[1 — a(T/To)™"(^''f'^'='') + ...] is in agreement with that obtained by Affleck and Ludwig [If. 
In particular, for the special case of = M, Agp = A^h = 1/2, we obtain a ^/T cor- 
rection as does the conformal approach [Q. Estimating the coefflcient a requires a knowl- 
edge of the temperature dependence of p^^ which is beyond the scope of the present pa- 
per. It is gratifying that for N = 2 the magnitude of the spin contribution to the rate, 
[(7rp)/(2r„^(0,0))] = 3tt^/[A{2 + M)^], (which is obtained after removing the 1/4 of the total 
rate due to potential scattering) agrees with the results of Affleck and Ludwig to within 
8% for all M > 2. 

Local Spin and Channel Flavor Dynamical Susceptibilities: The linear response to ex- 
ternal fields coupling to the impurity spin and channel quantum numbers can be easily 
calculated from the bubble diagrams for the spin and channel excitations, / and b, respec- 
tively. The leading and next-to-leading contributions to the absorptive part of the local 
spin and fiavor susceptibilities (per spin or channel degree of freedom), x'/^ = ImXsp/N and 
Xch = ^^Xch/M, are given by: 



Ludwig from conformal field theory 0]. 



Xch(^,0) 



Xspi^, 0) 




(11) 



(12) 
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Note that for N = M (which includes the important case N = M = 2) both suscep- 



tibihties reduce to the general form x"(u;,0) ~ y4sgn(ti;)[l — B\J\iu\/TQ + ...]. This lead- 
ing step function behavior in x" ^^id the associated logarithmic dependence in the real 
part, x'i^jT) ~ ~ln{max{uJ,T}/To), provides a possible connection with the marginal 
Fermi-liquid phenomenology of the high-Tc oxides |0| . This step function behavior was first 
noted in Ref. [4(b)]. For > M > 1 the real part of the spin susceptibility is constant 
(~ I/Tq) and the non- Fermi-liquid behavior is dominated by the divergence 
of the channel susceptibility. In the opposite limit, < M, the flavor susceptibility is 
constant (~ W^f^/To) and a |cI;|~(^ch-Asp) divergence occurs in the spin susceptibility. For 
N > M the system displays two-parameter universality: the channel fluctuations start dom- 
inating below a new energy scale, Tch ~ ToWj^^^^"'' '^ch)]^ with Tch << Tq foi N > M in the 
Kondo limit (Wch << 1). We note that Tch evolves into the "pathology" temperature below 
which non-Fermi liquid behavior ensues in the NCA treatment of the one-channel Kondo 
model |n],|12|. 



1/N Fluctuations: The functional integral formulation outlined above gives a natural 
framework for estimating the effects of fluctuations, and allows us to explain the remark- 
able fact that the saddlepoint exponents remain unchanged to all orders in The 
arguments are very much in the spirit of the argument usually made to justify the fact 
that perturbation theory gives the exact exponents characterizing the low energy behavior 
in conventional Fermi liquids. More precisely, all fluctuations can be incorporated 
into the renormalization of interaction vertices in all diagrams (/ and b self-energies, the 
single-particle Green function as well as all susceptibilities). We have checked explicitly to 
order that these vertex renormalizations modify the amplitudes but only give sublead- 
ing singular contributions. (For example, the sub-leading corrections to Sj and S{, behave 
as, lm6Ef{uj) ~ |u; — EqI^^"'' and Im5Sb(u;) ~ |u; — EqI^^^p, respectively.) In fact, the 
saddle point gives the exact low-energy singularities to all orders in as can be seen 

by considering arbitrary order diagrams in perturbation theory (around the saddle point). 
The appropriate propagators carry spectral functions which diverge as leu — i^ol^-^''"^ [see 
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Equations (6-9)]. As an illustration, consider a generic diagram contributing to Ilf, which 
contains L loops (thus L independent energy integrations), L fermion propagators, and 
L — 1 boson propagators. The most singular contribution behaves as \to — Eq]''''^^^ where 
Cb{L) = L + L{Af — 1) + (L — 1)(A{, — 1) = 1 — Afo since unitarity of the scattering ampli- 
tude requires Aj + = 1. This is indeed the behavior found at the saddle point. Similar 
arguments apply to Sj, the one particle Green's function and the spin and channel suscepti- 
bilities. It is natural to speculate that, in all systems with spin-charge separated, non- Fermi 
liquid ground states, the correct low energy behavior can be obtained on the basis of an 
appropriate (self-consistent) perturbation theory involving the "separated" spin and charge 
degrees of freedom. In analogy with conventional Fermi liquids, such perturbative argu- 
ments should be valid even when no obvious small parameters are available, provided no 
phase transition occurs to a Fermi liquid state through the "binding" of spin and charge. 

Crossover Effects: A related issue concerns the crossover to the Fermi liquid solution 
which becomes the correct ground state in the presence of a spin or channel symmetry 
breaking field, respectively Hgp or Hch- Our saddle point calculation leads to a crossover 
away from multi-channel behavior below respective scales, Tgp ~ Hlp^^^ ,Tch ~ H],^'^ . The 
corresponding crossover exponents, 0sp = 1 + 1/7, 'Pch = 1 + 7, are precisely those obtained 
from conformal field theory |^ . However, the Fermi liquid behavior which sets in well below 
the crossover scale Tgpi^ch) is outside the scope of the NCA. This can be traced back to the fact 
that the Kondo screening characteristic of the Fermi liquid fixed point involves the formation 
of a singlet bound state between a conduction electron and the local spin excitation, p!6| ; 
the residue of the bound state pole plays the role of the slave-boson mean field amplitude 
in the conventional large N approach to the Kondo problem. In the multichannel problem 
this Fermi liquid saddle point becomes possible only in the presence of a channel symmetry 
breaking field. 

Above we have shown that the limit — > oo, M oo with M/N = 7 fixed allows us to 
obtain the exact low energy behavior of the multichannel single impurity Kondo problem. It 
suggests the possibility of a phenomenological approach to the multichannel Kondo problem 
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and other impurity models with non-Fermi hquid ground states based on spin-flavor (or 
spin-charge) separated degrees of freedom. This can be regarded as the logical extension 
of Nozieres' classic discussion of the single channel Fermi liquid case to non-Fermi liquid 
situations. A detailed analysis of fluctuations about the multichannel saddle point, including 
the crossover to the Fermi liquid case will be described elsewhere. Moreover, this formulation 
leads to a natural extension to the lattice through the large-D treatment of correlated 
systems . 
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Figure Phase diagram in the N, M plane. The NCA is strictly controlled in the large 
M, N limit for M > 2, but gives the exact critical exponents for all N,M > 2. The 
universality class is the same for all lines of fixed slope 7 = M/N . 
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